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Abstract 

We present calculations of the stability of planar fronts in two mean field 
models of diffusion limited growth. The steady state solution for the front can 
exist for a continuous family of velocities, we show that the selected velocity 
is given by marginal stability theory. We find that naive mean field theory 
has no instability to transverse perturbations, while a threshold mean field 
theory has such a Mullins-Sekerka instability. These results place on firm 
theoretical ground the observed lack of the dendritic morphology in naive 
mean field theory and its presence in threshold models. The existence of a 
Mullins-Sekerka instability is related to the behavior of the mean field theories 

in the zero-undercooling limit. 
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I. INTRODUCTION 


Some of the most appealing and difficult patterns found in nature are formed by diffusion 
limited growth. A canonical example is solidihcation patterns [^,0, which, in addition to 
generating the delicate beauty of snowflakes, are important metallurgically as the source of 
alloy microstructure. Diffusion limited growth drives pattern formation in a wide variety 
of systems: examples include hngering patterns in the displacement of a viscous fluid [|^, 
electrodeposition patterns |^,Q, microbial aggregation |^,^], dielectric breakdown, flow in oil 
reservoirs, and possibly the wiring of neuronal dendritic arbors ||^. 

In recent years, much progress has been made on the problem of constructing a theory 
of pattern formation in dissipative systems P]. This theory explains the universal patterns 
seen in systems where a linear instability of a uniform ground state saturates at some hnite 
amplitude. In this context, diffusion limited growth patterns form an interesting example 
of a pattern forming system whose instability never saturates.^ 

The instability which drives pattern formation in diffusion limited growth is the Mullins- 
Sekerka instability. The physical nature of this instability can be seen by considering a solid 
seed particle immersed in a supercooled liquid bath. As the liquid solidihes around the seed, 
the latent heat released by the hrst order phase transition must diffuse away. An outward 
perturbation on the boundary of the growing solid can cool more easily, and so grows faster 
than neighboring regions, leading to an instability of the solid-liquid interface. Competition 
between the Mullins-Sekerka instability and other effects, such as surface tension, sets a 
basic length scale for pattern formation. As boundary features formed by the Mullins- 
Sekerka instability eventually grow to a size where they too are unstable, new structures 
grow on them, and highly complex fractal patterns are formed. 

Diffusion limited growth occurs in at least two broad morphologies: tip splitting growth, 

^Although situations can be created where some of the theory of weakly nonlinear pattern forma¬ 
tion can be applied, such as directional solidification |jl]. 
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as seen for example in Hele-Shaw flow in a radial geometry, and dendritic growth, as ob¬ 
served in solidification. In a radial geometry, tip splitting growth forms branching hngers 
whose ensemble averaged envelope is circular, while dendritic growth forms a concave en¬ 
velope lead by cusp-like tips propagating at some hxed velocity. The mechanism acting at 
the tips to select the actual behavior is understood: anisotropy in the surface tension (or 
phase transformation kinetics) acts as a singular perturbation at the tip to stabilize the 
dendrite, with the unique velocity being selected by a microscopic solvability condition [|l|. 
If anisotropy is absent, no propagating dendrite can be formed, and tip splitting growth is 
observed. | 

More recent work is on global aspects of these pattern forming systems. There is numer¬ 
ical evidence P-pT[| for a morphology transition between the concave asymptotic prohles of 
the dendritic growth morphology and the convex asymptotic prohles of the densely branched 
morphology (DBM). The work of Shochet et ah suggests that there is a discontinuity in the 
slope of the velocity vs. control parameter curve at the convex to concave transition, implying 
that the different morphologies are akin to thermodynamically different states. This kinetic 
phase transition hypothesis is bolstered by their suggestion that there may be coexistence 
between the branching and dendritic morphologies . 


One step towards theoretical understanding of this morphology transition and of dif¬ 
fusion growth patterns in general is to formulate an envelope model which displays the 
convex/concave transition of the underlying microscopic models. Such a model, based on 
the classic mean held models of Witten and Sander ||^ , has been proposed |T^-|17| . Since a 
model for the envelope cannot allow concavity (or the dendritic morphology) if it forms sta¬ 
ble planar fronts (by the kinetic Wulff construction discussed in Section B. it is necessary 
to look for a Mullins-Sekerka instability in the models discussed in |^, and to understand 
the nature of that instability. 


^Isotropic growth fingers can also be stabilized by growth in a sufficiently narrow channel Q. 
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In the remainder of this paper we will discuss several microscopic solidihcation models 
in general and review the behaviors of the mean held models (Section 0)- We will present 
some background on the geometrical construction of asymptotic shapes (Section B, the 
relationship between morphology transitions and the Mullins-Sekerka instability, and the 
velocity selection of planar fronts (Section [111 A|). Our calculation of the stability of planar 


fronts in mean held models will be presented in Section jlV 7 \| . The zero undercooling limit 
will be discussed in Section and our conclusions appear in Section 


II. MODELS AND MEAN FIELD THEORY 


Modelling ehorts for dihusion limited growth have taken a variety of forms. The direct 
continuum approach for a sharp interface is in terms of a moving boundary value problem 
1^. The held controlling growth [u, which is proportional to the undercooling — T) 
satishes the dihusion equation 

dfU = DV^u . ( 2 . 1 ) 


The boundary conditions are undercooling at inhnity and a melting point on the boundary 
set by a surface tension 7 (due to the Gibbs-Thompson condition) and phase transformation 
kinetics 


^ 1+00 

“Iinterface “ ~ (^'^n ■ 

The boundary motion is determined by the local gradient of the dihusing held 


( 2 . 2 ) 


W I interface ' ^ (2-3) 

where h is the normal unit vector of the interface. The displacement of a more viscous huid 
by a less viscous huid (the Sahman-Taylor problem) or how in porous media satishes similar 
equations, with the held u representing pressure in the huid |]^. 

More hexibility is ohered by phase held models, which specify a continuous transition 
from the solid to liquid phase Ill 8 |-p 0 l, allowing arbitrary topology of the interface and more 
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control over the boundary layer and the phase transformation at the expense of an additional 
held. 

Another approach is the diffusion limited aggregation (DLA) class of models ||T3| , pT| , 
which have generated a great deal of interest. This approach models diffusion limited growth 
by the discrete aggregation of random walkers, with the probability density of a random 
walk capturing diffusive behavior in a computationally efficient scheme. The random fractal 
nature of these aggregates motivate questions about the effects of noise. A phase held model 
with stochastic properties is studied in [^. Other models include the cell dynamical scheme 
^ and the dihusion transition scheme , both of which couple numerical solution of the 


dihusive held with some model for the phase transition in cells on a lattice. Noise is an 
inherent part of these models as well. 

The various diherent growth morphologies may all be observed within any particular 
model, depending on such parameters as boundary conditions, undercooling, surface tension, 
kinetic ehect, anisotropy, etc. The dendritic morphology occurs at hnite anisotropy, hnite 
undercooling; while single walker DLA corresponds to zero undercooling, zero surface tension 
and kinetic ehect. Tip splitting growth as seen in radial Hele-Shaw how and Sahman-Taylor 
hngers in a channel geometry occurs at zero anisotropy. 

A. Mean field theory 


Mean field models attempt to capture global morphology of dihusion limited growth. 
Naive mean field theory in open radial geometries at zero undercooling has a spherically 
symmetric solution which captures the mass scaling behavior of single walker DLA, giv¬ 
ing a mean field exponent of df = d—1, a. number believed to become exact as d ^ oo. 
However, this spherically symmetric solution is always unstable, leading to difficulties of in¬ 
terpretation. In addition, the naive mean field theory fails in one dimension or in a channel 
geometry, where fronts accelerate to infinite velocity in finite time . A phenomenological 
modification of the naive MFT has been proposed which fixes these problems and success- 
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fully explains several important morphological observations. One is that the envelope of 
ensemble averages of DLA in a channel geometry takes the shape of a Saffman-Taylor hnger 


and another is that anisotropic growth in a radial geometry can form the global den¬ 
dritic morphology |]^. There are now indications that the necessary threshold suggested 
phenomenologically in []H| can be generated by including the essential effects of noise in 
underlying model ||2^ . 

We begin our discussion of mean held theory by considering the Witten-Sander aggrega¬ 
tion model with hnite walker density. In this model, the diffusion held controlling growth is 
modelled by random walkers. When a walker becomes adjacent to the solid, it transforms 
into solid and advances the solid to that position. This sticking rule implies that the walker 
density u on the boundary is zero, giving a walker density adjacent to the boundary of aVu. 
This gives a growth velocity proportional to Vu. Comparison with (|2.2|) and (|2.3|) shows 
that DLA is a stochastic model for dihusion limited interface growth with no surface tension 
or interface kinetics. The random walker density will satisfy the dihusion equation, minus 
the walkers lost to the growing aggregate density p 


ii = DV^u — p . 


(2.4) 


The transition probability of an unoccupied site x is proportional to 

m(x, t)P(x, t) (2.5) 

where P(x) represent the probability of the neighboring sites to be occupied by the aggregate. 
We can write P(x) in the form 

+ ( 2 - 6 ) 

i 

where e* are lattice vectors to adjacent sites and i runs over the number of nearest neighbors. 
This formulation allows multiple walkers to occupy a single site, a feature not expected to 
ahect the universality class. 
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We derive a continuum mean field model by allowing u and p to take on any continuum 
values and identifying the sum of nearest neighbor terms in P as a lattice Laplacian. This 


results in the Witten-Sander mean held model for hnite walker density DLA: 

ii = DS/^u — p (2.7) 

p = Ku{p + a^V'^p) . (2.8) 

where the constant parameters P, K and a can be removed by setting x ax, t to?/D, 

and {p,u} — {p,u]D!a?K. In these units the equations become 

u = — p (2.9) 

P = m(p + VV) (2.10) 


with a boundary condition on the undercooling, or walker density at inhnity m|+oo = The 
total number of particles (aggregated and walking) p + w is conserved by the local dynamics. 

This mean held theory has no threshold for aggregate growth, so a small perturbation 
in the p held, in the presence of u, will grow exponentially. However, lattice simulations 
with discrete phase states, including the DLA model discussed above as well as others 
SHH, rave an implicit threshold for growth; growth at a site is disallowed unless a 


nearest neighbor site is fully occupied. Phase held models also have a threshold, in this 
case explicit in the structure of the form of the free energy, which specihes that the phases 
are at least metastable [^,|^]. This lack of a threshold in naive MFT leads to problems in 
channel or planar geometries, namely the inhnite propagation mentioned earlier. Inserting 
a threshold by hand, we have the phenomenological model introduced by Brener et al. ||14|| : 

ft = — p 

p = u{p^ + V^p) . 

Any function P(p) which vanishes faster than linearly as p —> 0 could be used in the place 
of p'*', for example P(p) = p0(p — A). 

The success of the 7 model motivates the question of how to obtain such a growth 
threshold in a more fundamental way. One possibility is to invoke the generic effects of 


( 2 . 11 ) 
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noise In microscopic theory, thermal noise at the tip drives sidebranching |]l[. In the 
Witten-Sander particle aggregation model, shot noise of particle diffusion and aggregation 
are a source of noise. Noise is known to shift the bifurcation point of nonlinear systems 
^, 111 . 


To sketch how noise may affect the form of the macroscopic mean held equations, recall 
the derivation of the naive mean held theory 


from the Witten-Sander rules. The 


transitions of particles on a nearest neighbor site from walker to aggregate are random and 
independent, and therefore governed by Poissonian shot noise. Including these huctuations 
by a Langevin noise source r] multiplicatively coupled with the held, we have 


p = Ku{p + a^W^p) + Ku{p + a?V^p) p 


( 2 . 12 ) 


Multiplicatively coupled noise can lead to stabilization of an unstable system, analogous in 
some ways to a parametrically driven pendulum [^. This stabilization is exactly what is 
required, as it establishes a lower threshold in p below which growth does not occur. A more 


complete discussion appears in |24 


III. THEORY OF KINETIC MORPHOLOGY 

In this section, we review some material concerning the growth of envelopes. A system 
with a transition layer that does not spread will in the long time limit be described by a 
sharp interface model. We parameterize the interface between the two phases by a curve 
x(cT, t) and write the evolution of this curve in the form 

dx , , 

— = v{a,t,[x]) (3.1) 

where v can in general be a complicated functional depending nonlocally on the entire curve 
X. We are interested in envelopes which form a well dehned shape in the asymptotic time 
limit, that is assume a scaling form 


T(cr, t) = R(t)r{a) . 


(3.2) 



















Since velocity parallel to the interface corresponds merely to a redefinition of the parame¬ 
terization a, we take v = Vnfi. From ( |3.1|) and (|3^) we get 

fix 

— = R{t)f{a) = Vn{(J, t)h (3.3) 

If the envelope scales, the time dependence will factor out of Vm leaving us with 

Vn{(^) = r{a) ■ n . (3.4) 


One may imagine several possibilities for the nature of For example, in a sharp interface 
model for diffusion limited growth, Vn will depend on the whole structure and history of 
the envelope in some nonlocal way. Alternatively, it may be that v depends only on purely 
local properties of the interface, such as the normal direction, curvature and so on. In the 
asymptotic limit, the envelope locally appears to be flat, leaving the normal direction as the 
only remaining variable that v can depend on. In this case, one can solve for the envelope 
shape explicitly given Vn{n) for all angles of flat front propagation. This is done by the 
kinetic Wulff construction to be discussed next. 

The equilibrium shape of a crystal has a solution in the form of the Wulff construction 
which dates back to the turn of the century p7| , p8| . Given the anisotropic form of the surface 
energy 7 ( 6 *) for each crystal direction 6, minimizing the total surface energy (at hxed area) 
leads to the following construction: for every normal direction, draw a line perpendicular 
to the normal at a distance from the origin proportional to surface energy of that direction. 
The inner convex hull of these lines is the equilibrium crystal shape. The resulting crystal 
may be faceted, smooth (if it is above the roughening transition) or may have both facets 
and smoothly curved surfaces. 

If a crystal is grown under such conditions that the normal velocity depends only on the 
direction of the surface normal (for example, in a reaction limited regime), then the steady 
state shape is again given by the Wulff construction, with normal velocity substituted in for 
surface energy. It is clear that such a steady state cannot be concave [^,0. Because the 
envelope of the dendritic morphology is concave |n|, a geometrical model cannot capture 
the dendritic morphology. 
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Although the kinetic Wulff construction is not applicable to the nonlocal diffusion limited 
regime, it is important to understand its consequences for general modelling of the morphol¬ 
ogy transition. A model which forms stable planar fronts at any angle in two dimensions 
cannot capture the dendritic morphology. In the asymptotic limit, a scaling form will arise, 
composed of approximately planar fronts. The velocities of these planar fronts as a function 
of angle can be found by marginal stability. The asymptotic shape will therefore be given by 
the kinetic Wulff construction, and will be convex. So a mean field theory which always has 
a stable planar front cannot model the dendritic (concave) morphology, and cannot describe 
the morphology transition. In fact, in the convex regime, the planar front has to have an 
instability analogous to Mullins-Sekerka instabilty. It is for this reason that we have tested 
for the presence of a Mullins-Sekerka instability in mean field models. 

In our calculations we have not included a term such as 




(3.5) 


dy^ 

to explicitly include anisotropy. Although such an anisotropy term is of course a necessary 
ingredient for dendritic growth, we have neglected it for simplicity. Its effect will be small; its 
function at the level of our linear stability analysis would merely be to pick out a preferred 
direction for the instability to grow. This point of view is verified in 2d simulations, which 
show explicitly the existence of dendritic growth in exactly those systems with a Mullins- 
Sekerka instability in isotropic planar fronts [j^ . 


A. Marginal stability 

To find the front velocity v{9) used above, it is necessary to find the velocities of planar 
fronts in the desired model. The mean field diffusion growth equations derived earlier have 
the generic property (shared with other nonlinear diffusion equations) that there exist steady 
state solutions of any velocity. This property is connected with the unstable nature of the 
invaded state, so it does not apply to models with a strict cutoff, although it does apply to 


10 







the 7 model with its power law cutoff. The hrst task is therefore the solution of a velocity 
selection problem. Marginal stability theory states that, for typical initial conditions, the 
selected velocity is the lowest velocity for which localized perturbations are stable in the 


moving reference frame |31|J^ . Here, “localized” means perturbations whose asymptotic tails 
decay no slower than the front itself, and are therefore smaller than the front everywhere. 
Marginal stability is on a hrm foundation for nonlinear diffusion problems with one variable. 
Marginal stability falls into two classes: linear marginal stability, where computing the 
stability of the asymptotic tail is sufficient, and nonlinear marginal stability, where the 
(linear) stability of the entire front needs to be considered. In one variable velocity selection 
problems, such as the Fisher equation u = u{l — u) + V^u, a parameter counting argument 

i.e. 


indicates a continuum of unstable modes in the regime where the front oscillates |31 


where c < 1/2. In problems involving nonlinear marginal stability, it is the appearance of a 
discrete unstable mode which sets the marginally stable velocity. 

There is a qualitative distinction between fronts selected by linear and nonlinear marginal 
stability. In either case the front is formed as a balance between the growth of the density 
and its diffusion which catalyzes further growth. In the linear case, the growth of the front 
at its leading edge is sufficient to catalyze its continued propagating, in this case we speak of 
the front being “pulled” by its leading edge. In the nonlinear case, the linear growth rate of 
the leading edge is insufficient to maintain the propagation velocity, and the larger nonlinear 
growth in the body of the front is necessary to catalyze the continued propagation. In this 
case, we speak of the front being “pushed” by the growth behind the front . 


To compare with the one dimensional marginal stability theory we compute, for the naive 
(7 = 1) model, the asymptotic behavior of the front at z = +cx). This can be found by 
setting M = A in the equations for a steady state front (see Section giving 

- cp'= A{p + p") . (3.6) 


This has solutions p ~ exp \±z where \± is given by 




4A2 


1 . 


(3.7) 
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For c < 2A, these modes oscillate, which can be prohibited on physical gronnds by the 
interpretation of p as a density. The parameter connting argnment analogons to the nonlinear 
diffnsion case is complicated by an additional variable and an additional free parameter (A) 
bnt the resnlt is the same: namely, that we have a continnnm of instabilities in the regime 
where the helds acqnire oscillating tails, predicting a linear marginally stable velocity of 
c = 2A. 

For 7 > 1, we expect similar behavior — the selected velocity will be the slowest stable 
velocity, which will again be the lowest velocity for which p does not cross zero. In this case, 
however, this will not be given by the behavior of the asymptotic tail. The linear growth 
rate for 7 > 1 is zero, and therefore any front formed mnst be pushed, with a velocity given 
by nonlinear marginal stability. 


IV. PERTURBATION OF MOVING FRONTS 


We will write out the equations describing perturbations of fronts satisfying (| 2 . 11 | ). We 
consider a planar geometry, so that the only spatial dependence is x dependence, and we 
will add transversal dependence back in later. Working in a moving frame z = x — ct and 
looking for a steady state solution = 0 we have 


— cu! — u" — CfJ 

(4.1) 

- c,l = «(F(p) + p") 

(4.2) 


where primes indicate z derivatives. By integrating and using the boundary conditions 
p = A, u = 0a.tz = —00 the hrst equation becomes 

u' = c(A — p — u) . (4.3) 


By integrating forward from z = —cx), a steady state front can be found for any desired 
velocity. 

The base solutions may be found by numerically integrating the equations. Initial condi¬ 
tions at the left are chosen to assure the boundary conditions u —0 and p ^ A as z ^ — 00 . 


12 





Since the p field may become negative, we use sgn(p)|p|'>' as our threshold function. Next 
we do perturbation theory on our equations in the moving frame 


u — cdzU= — p + cdzP (4.4) 

p - cdzp= u{F{p) + VV) (4.5) 

where we include transversal dependence in our perturbation 

u= vP{z) + {z) (4.6) 

p=pO(;,) + e‘^W(^) (4.7) 


and look for eigenfunctions and p^ which go to zero as > ±cx). Linearizing, we get 


1 1 / 721 1 1 

UU — CU = —K U -\- U — up + Cp 


(4,8) 


u)fi' - cp'' = u'’(c)(F'(p°(c))p' - k^p' + p'") 

W(F(p"(z)) + p“'P)) (4.9) 

Since there is little hope of solving this eigenvalue problem exactly, we discretize it in 
a box and diagonalize it numerically. To put this in a form which is easy to discretize, we 
write it as 


D^ + cD- (F(pO) + pO" + e) 

-u^D^ - u^{F'{p^) - e) 



= u 


(F(p«) + p»”) 

u^D^ + cD + u^{F'{p^)-e) 





where D represents the finite difference approximation of a derivative expressed as a tridi¬ 
agonal matrix. 

Now, we also want our perturbation to decay as fast or faster than the base solution 
a,t z = -|-cx). To do this, we can rewrite our eigenvalue problem in terms of new variables, 
having divided out the asymptotic fall-off of the base solution tails on the right, after which 
a zero-slope boundary condition on the right enforces the localized perturbation condition. 
Explicitly, if we let = fu{z)ul^g^, and ph^ = fp{z)Pn^^, where the / functions have the 
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same dependence at large z as the base solution and approach a constant at the left, then 
we have 


Dv} = DUz)ul„ = Uz)[§- + D)ul^ (4.11) 

Ju 

and 

Uz)S + Dful^ 

Ju 

=/.W( j + . (4.12) 

jtf pf pff 

Making the replacements D —> D + y and D'^ —> ^ in the perturbation 

equations (f4.10|) , we have 

D^ + (c + 2f )D + (-F(pO) - pO" _ f + cf) 

^ Ju ' \ \i / I / 

+ rn 

- 2^u’>D - u^fF'fp") - e + !f)} 

Ju Jp \ \I / Jp ' 

+ (c + 2^u^)D + u^Fjp^) - P + ^) + 

The functions / are computed from the base solution for 7 > 1 and chosen analytically 
for 7 = 1. The growth rates of possible perturbations are just the eigenvalues of this 
nonsymmetric matrix, which can be found numerically by standard techniques. 



1 


1 


u 


u 



= U! 







p 


p 


(4.13) 


A. Numerical results 

We hrst display the results relating to velocity selection by marginal stability. The 
behavior of the largest growth rate a;max near Cms = 2 A for 7 = 1 is shown in Figure |I]. As 
expected, we see that the growth rate for the most unstable mode drops to zero as c passes 
through 2A, implying a selected velocity of 2A. This velocity selection has been verihed 
in direct numerical simulations. At 7 = 1, then, the front is pulled, and linear marginal 
stability holds. 

The prohles of the u and p helds at a velocity smaller than Cms for 7 > 1 are shown in 
Figure |^. At this velocity, the p held crosses zero to become negative. This creates a well 
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allowing a discrete unstable mode, analogous to a bound state in a well in quantum wave 
mechanics. For 7 > 1, we again expect the velocity to be given by the velocity of the slowest 
stable front. This is now a global stability calculation for the velocity of a pushed front. 
This velocity selection is verihed in Figure |^, which compares the velocity calculated from 
numerical diagonalization at 7 = 1.2 with the velocity as measured in simulations. It is clear 
that nonlinear marginal stability theory is predicting the correct velocities as measured in 
simulations, and that the 7 model gives pushed fronts. 

The nature of the instability for 7 > 1 is indicated by the spectrum shown in Figure 
As c passes below Cms, a single discrete mode moves above zero. This is in contrast to the 
7 = 1 case, where the top of a continuum moves above zero, but it is natural in light of 
the fact that p decays as a power law and cannot have the oscillating tail that generates 
the continuum in the 7 = 1 case. Figure ^ shows a quadratic extrapolation to zero velocity, 
which occurs for 7 = 1.2 at A = A* ~ 0.055. This dehnes A*, the lowest undercooling 
at which the model can form a propagating front. Nonzero A* is the essential feature of 
threshold mean held theories, and is not found in naive mean held theory, for which (since 
Cms = 2A) A* = 0. The transverse stability calculation is done for the selected planar front, 
and this calculation gives the required front velocities. 

Our main results concern the existence of a Mullins-Sekerka instability in these 
models. Again we restrict ourselves to perturbations which decay at least as fast as the 
front. Figure ^ shows the growth rate uj of the most unstable mode in naive mean held 
model plotted as a function of the transverse wavenumber k for a number of undercoolings 
A. There is no sign of any hnite wavenumber instability. 

Figure ^ shows the largest growth rate 00 vs. k for a number of A’s in the 7 model at 
7 = 1 . 2 . There is a hnite k instability for A below A** ~ 0.29. Figure shows a typical 
unstable spectrum. The characteristic length scale of the instability is comparable to the 
width of the front. Comparing with Figure note that there is no sign of any velocity slope 
discontinuity at the point A** ~ 0.29 where the Mullins-Sekerka instability disappears. This 
argues against the morphology transition being some kind of nonequilibrium phase transi- 
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tion. These results demonstrate explicitly what was already indicated in a threshold 
is necessary to generate a Mullins-Sekerka instability and allow concave dendritic envelopes 
to be formed, and that the 7 threshold model in particular both has concave envelopes and 
a Mullins-Sekerka instability. 


V. THE LIMIT OF ZERO UNDERCOOLING 


Undercooling is an important control parameter in diffusion limited growth in general 
and the only control parameter in many walker DLA. The diffusion transition scheme ||10|| , 


for example, becomes equivalent to the Ising model when the undercooling is maximal. 
Many walker DLA is equivalent to the Eden model when the density of walkers equals one. 
The single walker DLA model of Witten and Sander |T^ corresponds to the opposite limit 
of zero undercooling. As undercooling is changed, morphology transitions can be observed. 
Varying the undercooling from high to low in the cell dynamical scheme, Liu and Goldenfeld 
m found transitions from dendritic growth aligned with crystalline axes to a more or less 
isotropic state (interpreted as the dense branching morphology), and then back to dendritic 
growth, this time with dendrites aligned between the crystalline axes. 

Naive mean held theory for single walker DLA takes the following form: 


V^u - p- u{p + V^p) , 


(5.1) 


where m|+oo ~ Jf, i-e. there is some hxed hux of u at inhnity. This is the form in which 
DLA mean held theory hrst appeared [0, however it is not a simple limit of the hnite A 
theory. To see the connection, we start with (|2.7|J2.8|) and put u ^ u/A, t ^ At, obtaining 


p = u{p + VV) 
Ait = — p . 


(5.2) 

(5.3) 


The argument is now that as A —>■ 0, we have = p, which is the mean held model at 
zero undercooling quoted above. The question is what happens to the boundary conditions. 
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We can relate the fixed fiux boundary conditions of the usual zero undercooling model to 
the A —> 0 limit of the finite undercooling model by considering the system in some box 
with a finite length L. After the rescaling above, the boundary condition on u is m|+oo = 1- 
In a system of length L, this corresponds to a supply of fiux J ~ 1/L. Taking the system 
size to infinity shows that the direct A —0 limit of the finite undercooling model is the 
zero undercooling model with zero fiux. This somewhat odd condition is due to the nature 
of the Laplacian model: it is modelling infinitely slow growth, so it must introduce a new 
artificial time scale (given by the boundary fiux J) on which growth formally proceeds. 

At first sight, the lack of a Mullins-Sekerka instability in naive mean field theory seems 
surprising, as it seems to contain all the ingredients which lead to a Mullins-Sekerka insta¬ 
bility in other diffusion limited growth systems. In fact, there is a connection between the 
existence of a Mullins-Sekerka instability and the behavior of fronts in the zero undercooling 
limit alluded to earlier. Specifically, any system with propagating fronts in the zero under¬ 
cooling model should, in the finite undercooling version, 1 ) have a lowest undercooling A*, 
below which propagating fronts cannot form and the boundary region spreads forever and 
2) have some region of instability above A* and below some value A** where the system 
can form propagating fronts in one dimension but where those fronts exhibit a transverse 
Mullins-Sekerka instability. These effects are both due to the existence of a threshold for 
growth in the model. 

For a zero undercooling model (such as the 7 model) which possesses a travelling front 
solution, there is some selected value p° for the density of the aggregate left behind by 
the front [Q. This value is independent of the boundary fiux J, which can be scaled 
out by rescaling time. Models such as the naive model do not have such an aggregate 
density; as shown in via a scaling solution the frozen cluster scales with distance as 
1/x^, corresponding to a “fractal dimension” oi df = —1. The constant fiux of walkers into 
the system forces the front to advance at an ever-increasing rate, reaching infinite velocity 
in finite time. 

Consider now a finite undercooling model as we approach the limit of zero undercooling. 
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The aggregate density must equal A by conservation of matter. In the zero undercooling 
version, we know that if a steady state front forms it must leave behind an aggregate density 
of some hnite p°. As A is lowered below p° in the hnite undercooling version, it will no longer 
be able to form a travelling front, because its asymptotic aggregate density is insufficient to 
support a travelling front in the zero undercooling version. Therefore a minimum undercool¬ 
ing A* for travelling fronts in the finite undercooling model is a signal of travelling fronts 
with an asymptotic walker density p° in the zero undercooling version, and vice versa. We 
have compared A* and p° as measured in simulations for various values of 7 , and found a 
rough correspondence consistent with this argument and the quality of our numerics. 

Next we turn to the question of the existence of a Mullins-Sekerka instability. The 
existence of a lower threshold for propagation A* indicates that, in some sense, the front 
has a “diffusion problem” at low undercooling: it has difficulty attracting enough walkers 
to grow. In the regime just above A*, we should not be surprised that the system exhibits a 
Mullins-Sekerka instability, as an outward perturbation can grow by shielding neighboring 
regions. The region below A* in which the system cannot form a propagating front is a 
region in which there must certainly be a Mullins-Sekerka instability (compare with the 
instability of a motionless interface P3| , |2(][| ). This instability persists in the moving front 
as well, at least up until some velocity fast enough to damp it; a phenomenon well known 
in the standard models of crystal growth. In a model with A* = 0, there is no diffusion 
problem, and no effective shielding can take place. We found stability in this case. 

VI. CONCLUSION 

The diffusion limited growth problem is a complex and widely studied exemplar for 
pattern formation far from equilibrium. The patterns formed can be quite regular, such 
as the organized crystals of the Dendritic morphology, or surprisingly complex, such as the 
fractal structures of isotropic single walker DLA. 

To understand the stochastic behavior of DLA like patterns, statistical methods have to 
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be used. For other microscopic growth models, like the Eden model, ballistic deposition, 
molecular beam epitaxy etc., this statistical approach in the form of nonlinear Langevin 
equations is quite successful in determining the scaling properties of the pattern and the 
growth process. Notable among them is the success of the noisy Burger’s equation (KPZ 
equation) for studying kinetic roughening processes such as the Eden model. While these 
models are all local growth models, DLA is non-local, which makes the construction of a 
proper held theory for DLA much harder than the local models. 

Mean held models, as a hrst step towards a fully renormalized theory, play a central role 
in the equilibrium theory of critical phenomena, and the corresponding theory in nonequi¬ 
librium systems is becoming better developed. In spatio-temporal chaos, for example 
in the Kuramoto-Sivashinsky equation , and in dihusion-limited reaction systems 


statistical averaging and mean held modelling are being used to obtain nontrivial macro¬ 
scopic behavior from microscopic models. Work on mean held models of dihusion limited 
growth thus hts into a larger context of ehorts to understand pattern forming systems far 
from equilibrium, and must be understood before steps going beyond mean held theory can 
be taken. 

Unlike most of the other cases, where the behavior of the mean held theory (the deter¬ 
ministic part of Langevin held theory) is rather straightforward, the mean held theory of 
dihusion limited growth is more involved. As we have shown before, a proper cutoh has to 
be introduced to avoid singular behavior of the mean held theory in low dimensions. Our 
modihed threshold mean held theories can exhibit the morphology transition, since they 
have the appropriate instability for the dendritic morphology. 

In this paper, we have made a precise stability calculation for the threshold mean held 
theory and for the hrst time established the existence of a Mullins-Sekerka instability for 
models with 7 > 1. This instability puts the widely observed transition between the den¬ 
dritic morphology (with concave envelope shape) and the dense branching morphology (with 
convex envelope shape) on hrm theoretical ground. We have also demonstrated the absence 
of such instability for the naive mean held theory without a cutoh. The drastically diherent 
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behaviors between the model with and without a cutoff can be related to their different 
behavior at zero undercooling and strongly reconhrms the existence of some sort of cutoff 
for the correct mean held theory. 
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FIGURES 



FIG. 1. Marginal stability plot for naive (7 = 1.0) mean field theory showing growth rate of 
most unstable mode vs. c. A = 0.25. Below c = 2A = 0.50, the p field crosses zero and has an 
oscillating tail.The displacement of the zero mode away from zero is due to the nature of the front 
tails right at c = 2 A, and does not occur for 7 > 1 . 



FIG. 2. p and u fields at a velocity below the marginally stable velocity. Note how p crosses 
zero. A = 0.25, 7 = 1.2, c = 0.15. At this 7 and A, the selected velocity Cms = 0.216. 
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FIG. 3. Front velocity as a function of A, as calculated by marginal stability and measured by 
simulations. The marginally stable velocity is the smallest velocity for which the front is stable in 
its own frame against local perturbations. 
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FIG. 4. Portion of the spectrum below the marginally stable velocity. 7 = 1.2, A = 0.25. Two 
sizes of grids are shown, to give some indication of the nature of the modes in the continuum limit. 
Most modes not shown lie on or close to the negative real axis. 



FIG. 5. Growth rate of most unstable mode as a function of transverse wavenumber of pertur¬ 
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FIG. 7. A portion of the spectrum for 7 = 1.2, A = 0.20, and k = 0.10. 
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